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Abstract. The results obtained in paper I are used to study possible relationships between the truncation radius of stellar discs 
in the NIR and structural parameters of the galaxies. The NIR truncation radius is larger for brighter galaxies, being proportional 
to W c m with c « 3/2, and with V m being the asymptotic rotation velocity at large radii (when the rotation curve becomes flat), and 
is lower for higher wavelengths. When it is normalized to the scalelength, the truncation is an increasing function of the central 
surface brightness and is lower for late type galaxies, although these correlations are weaker. These relations are in agreement 
with the scenario of magnetically driven truncations. 



Key words, galaxies: structure - infrared, magnetic fields 



1. Introduction 

Paper I (Florido et al. 2005) presented the result of the ob- 
servation of 18 edge-on spiral galaxies in the NIR in order to 
study the truncation of the stellar disc, the main objective be- 
ing to describe the truncation curve. This curve is defined as 
the progressive decrease of the photometric profile taking place 
at large radii with respect to the mean exponential profile ob- 
tained at moderate radii. In this paper we study the relationship 
of truncation to other parameters of the galaxy. We investigate 
some key properties that could restrict some of the theoreti- 
cal scenarios that have been proposed: what types of galaxy 
show truncation and what relation exists between the trunca- 
tion, radius, mass, luminosity, rotation velocity at large radii, 
type, central surface brightness and age of stellar population. 
The results obtained from this study will be merged with those 
from a previous one (Florido et al. 2001) in order to obtain a 
larger sample. No merger with other observations in the optical 
has been attempted, as the assumption of the same behaviour 
for optical and NIR truncations is not theoretically justified. 

In early works by van der Kruit (1979) truncations were 
found to have sharp cut-offs. In paper I and in Florido et al. 
(2001) truncations were found to be smooth. In these latter two 
papers it was shown that the smooth behaviour of truncations 
was not due to projection effects, but was real, existing even 
in deprojected profiles. In some recent papers (e.g. Pohlen et 
al. 2002, Perez 2004) this description is replaced by a two- 
exponential profile with a sharp elbow. There is increasing ev- 
idence favouring the two-slope description, at least in the opti- 
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cal. We are not interested here in the shape of the profile (see 
paper I); it may even be that the differences between authors 
are due to a true difference between the NIR and the optical. 
For the present, it is not our goal to determine whether this dis- 
crepancy is real or is a result of the low signal/noise at these 
very large radii, especially important in the NIR, and we refer 
to a "truncation (break)" or simply "truncation". 

In paper I it was proposed that most galaxies in the NIR 
seem to have a truncation curve that is proportional to (R - 
R tr Y n where n « 1, R is the radius and R, r is the truncation ra- 
dius at which the truncation seems to be complete, as obtained 
by extrapolation of the data. 

2. Theoretical models 

There are at present at least four alternatives to explain trunca- 
tions (or breaks): a) the collapse model, b) the threshold model, 
c) the magnetic model and d) the interaction model. 

According to the collapse model, which was introduced 
by van der Kruit (1987), truncations were formed in the col- 
lapse of the protogalaxy, and would correspond to the maxi- 
mum angular momentum of the protogalaxy. The gas beyond 
the truncation was accreted in later phases of galactic evolu- 
tion. Ferguson and Clark (2001) presented a model in which 
the truncation/break is also a result of the initial conditions in 
the gas. 

The threshold model was introduced by Kennicutt (1989). 
Here, when the gas density falls below a certain threshold, 
star formation should cease. A recent model by Schaye (2004) 
improved the dynamic mechanism by which the star forma- 
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tion threshold is established. Van den Bosch (2001) presented 
a hybrid model combining the collapse and threshold models. 
Elmegreen and Hunter (2006) reproduce two-slope profiles in 
a threshold model which includes several local effects of SF. 

Kregel and van der Kruit (2004) have shown that the thresh- 
old model is better supported by the observations than is the 
collapse model. The model by Schaye (2004) predicts an in- 
crease of R tr lh for small radial scale lengths and an increase of 
Ri r /h for higher central intensities and for low surface bright- 
ness discs. 

For an extensive discussion of the different models, see 
the review by Pohlen et al. (2004). See also the discussion in 
Battaner et al. (2002). 

In the magnetic model of truncations (Battaner et al. 2002), 
a centripetal magnetic force acting on the gas makes an im- 
portant contribution to the redistribution of stars once they are 
born. If the gas rotates under the equilibrium of gravity plus 
magnetic forces versus the centrifugal force, and if the mag- 
netic force is suddenly suppressed in the star formation process, 
the newborn stars would migrate and eventually escape. Some 
important modification of the mass distribution of the stellar 
density is therefore unavoidably expected at the region where 
the magnetic force becomes non negligible in the dynamics of 
the gas. 

The interaction model is only rarely defended (Noguchi & 
Ishibashi 1986). Schwarzkopf & Dettmar (2000) found similar 
properties of truncations for isolated galaxies and for interact- 
ing galaxies. However, close inspection of the large amount 
of data provided by these authors reveals some differences. 
Interacting galaxies have a ratio R, r (R)/R, r (K) that is always 
larger than one (R is red, K in the NIR), i.e. the truncation ra- 
dius is larger at shorter wavelengths (with a mean ratio of about 
1.6); however, in non-interacting galaxies, this ratio is closer to 
unity. Thus interaction effects should not be ruled out. 

3. Results 

Paper I presented new observations of edge-on galaxies ob- 
served in the NIR. In all, 18 galaxies were observed, but 5 of 
them were too noisy, 4 untruncated (within observational lim- 
its) and 2 uncertain. Therefore, only 7 galaxies were available 
to study truncations. Five of them, namely NGC 522 (Sbc), 
MGC-01-05-047 (Sc), NGC 2862 (SBbc), NGC 3279 (Scd) 
and NGC 5981 (Sc), were very clearly examples of this phe- 
nomenon. There are also two galaxies in which the truncation 
is not evident but that can be included in our study. These are 
NGC 781 (Sab) and NGC 3501 (Sc). We also use three trun- 
cated galaxies from Florido et al. (2001):NGC 4013 (Sb),NGC 
4217 (Sb) and NGC 6504 (Sbc). Therefore the subsample of 
truncated galaxies adopted in this study is formed by 10 galax- 
ies, 7 from paper I and 3 from Florido et al. (2001). As we have 
data for both sides of the galaxies, and measure at least in two 
filters (J and K s ), we have 2x2x10 = 40 points for our statistical 
subsample. 

However, these 40 data points are not independent. Because 
of the azimuthal symmetry of the galaxies, the results obtained 
at both sides should be similar (as, in fact, they are) and there- 
fore correlated. Moreover, the closeness of the filter wave- 
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Fig. 1. Ratio of truncation radius to scalelength (R rr /h) 
versus morphological type T (from the LEDA database, 
http://leda.univ-lyonl.fn. Symbols are: open square for K s and 
filled square for J. In the first panel we plot R tr /h when R lr 
and h are obtained by the same filter. In the second one we plot 
R, r /h(K s ) using h(K s ) for this ratio with each filter. 

lengths, for J and K s (for a few galaxies, also H), does not 
permit us to consider the J and K, points as really indepen- 
dent, except when considering a differential colour behaviour. 
The fact that we measure in two filters and on both sides per- 
mits us to work with smaller errors, but strictly speaking we 
have only 10 points for our sample, i.e. the number of galaxies 
selected. Despite the smallness of the sample, we present the 
more complete data for truncations in the NIR. 
A) Truncation versus Hubble type. 

At first glance, even for such a small number of galaxies, it 
seems that truncations are found more often in late type galax- 
ies. 

This is a suggestion and is not based on statistics, as the 
number of galaxies is too low and many more observations 
would be required. However, we take note of this potential dis- 
crimination, as it could be coherent with other findings shown 
later. This suggestion should be considered with caution, also, 
because the selection of our samples was not based on physical 
criteria, but rather considering appropriate angular size, large 
galactic latitude and date of observations. Therefore the possi- 
bility that a bias may be responsible cannot be ruled out. 

We looked for a relation of R tr (J)/h(J) versus T, where R, r 
is the truncation radius, h is the radial scale length and T is the 
galaxy type. J and K s are the filters used (Fig.l). The following 
tendency was observed: the larger the value of parameter T, 
the smaller that of R, r (J)/h(J) is. This fact is consistent with 
the above suggestion that truncation is a phenomenon mainly 
affecting late-type spirals. 

For some relationships, if we want to adopt /iasa length 
representative of the size of the galaxy, then h(K s ) better rep- 
resents the "true" radial scale length, as it is less affected 
by extinction and better represents the old stellar population. 
Therefore, we have plotted R, r (J)/h(K s ) and R„(K s )/h(K s ). In 
this case, the relationship seems to be reinforced. 

Therefore, truncation could be a phenomenon that mainly 
affects late type spirals. This relation is also found when the 
large sample of Schwarzkopf & Dettmar (2000) is used. This 
suggestion, although attractive, remains to be confirmed. 
B) Colour dependency of truncations 
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Table 1 reveals an evident property: R tr (J) > R, r (K s ). This 
is so even if the differences are small, because in all galaxies, on 
all sides, we find R tr (J) > R tr (K s ), with the exception of NGC 
6504. Therefore, there is a dependence of R tr with colour. This 
number of values is not small but we have applied three tests 
in order to discard an effect of chance: the T test, the Wilcoxon 
test and the sign test. In all three tests, the null hypothesis of 
null difference between R, r ( J) and R tr (K s ) was rejected. Taking, 
as usual, a contrast size of a — 0.05, in every case the signif- 
icance level (P-value) was much less than 0.05, being 0.004 
for the T test. Therefore, R tr (J) is unambiguously higher than 
R tr (K s ). This significant difference informs us about a relation 
with the age of the population. 

With respect to relative values, taking the radial scale- 
length as the unit, we obtain < (R tr /h) > j- 3.12 with <x = 0.85 
and < (R tr /h) >k s = 3.21 with cr = 0.93. As mentioned above, 
h(K s ) better represents the size of the galaxy. Then, we ob- 
tain < R lr (J)/h(K s ) >= 3.57 with cr = 1.06. Of course < 
R tr (J)/h(K s ) > is significantly higher than < R, r (K s )/h(K s ) >. 

These values are in reasonable agreement with other val- 
ues found in the literature (Pohlen et al. 2002, Barteldrees & 
Dettmar 1994). For comparison, it should be taken into account 
that some authors give Rb r or the point where the break takes 
place and others where it ends. Furthermore, our data in the 
NIR do not have to coincide with optical data, the differences 
being important. 

C) Truncation and rotation and luminosity. 

Figure 2 clearly shows a relation between R tr and V m , the 
asymptotic rotation velocity (taken from LED A). R tr (measured 
in kpc) is an increasing function of V m , which is emphasized in 
the log-log plot, showing a relation of the type R, r oc V c m . A re- 
lationship between R lr and V m was first obtained observation- 
ally by Pohlen et al. (2004). This is clearly supported by our 
present findings, in another wavelength range, with less scat- 
tering and with the mathematical function being determined 
quantitatively. This fact can be used to restrict the number of 
theoretical models, as we will see in the discussion. 

As different truncation values were detected for J and for 
K s , we have analyzed this relation (log/? (r , log V m ) for the two 
filters independently (Fig. 2). Standard Pearson correlations are 
used because the relation appears to be fairly linear in log space 
with classical significance tests. For [log R tr (J), log V m ] the cor- 
relation coefficient is 0.83, the slope is 1 .54 with error 0.28 and 
the log Rt r at the origin is -2.45. The P-value is 0.001, much less 
than 0.05, and therefore the null hypothesis of the absence of a 
linear relation between both quantities is clearly rejected. For 
[log R, r (K s ), logV m ] the correlation coefficient is again 0.83, 
the slope is 1.58 with a typical error of 0.29 and the value of R, r 
at the origin is -2.61. The P-value is 0.002, again much lower 
than 0.05. For both filters, the slope of the regression line (1.54 
+ 0.28 and 1.58 + 0.29) can be approximated to the value 1.5, 
which is more useful when looking for simple low-power rela- 
tions. 

A relationship physically connected to this one (R tr , V m ) is 
the relation (R, r , L) with L the luminosity of the galaxy. This 
is shown in Fig. 3, where we plot \ogR tr versus the absolute 
B -magnitude, showing a clear relation, with R tr being larger 
for brighter galaxies. A statistical analysis is unnecessary as 
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Fig. 2. Log of truncation radius versus log of maximum veloc- 
ity rotation. Symbols as in Fig. 1 
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Fig. 3. Truncation radius versus absolute B-magnitude. 
Symbols as in Fig. 1 . 



this relation is a consequence of the previous one, taking into 
account the Tully-Fisher relation. 
D) Other relationships. 

Fig. 4 plots R, r versus h and shows a clear correlation be- 
tween the two parameters. Thus, more extended galaxies have 
larger truncation radii, in agreement with Kregel (2003). Larger 
galaxies have lower truncation radii when measured with the 
radial scale length as the unit. This relation seems to be non- 
linear. 

Figure 5 shows R, r /h versus /j. (deprojected central surface 
brightness). This relation has been studied by Kregel & van der 
Kruit (2004), as mentioned above. Although the scatter of the 
points is much higher, we see that R tr /h seems to be a decreas- 
ing function of /v The correlation coefficient for [R tr (J), fi ] 
is only 0.39 and is the same for [R, r (K s ), fi a ]. Therefore, the 
non existence of this relation cannot be statistically rejected, 
as the P-value is (in both cases) 0.091, larger than the usual 
value of 0.05. The correlation coefficient for [R, r /h(K s ), /i ] is 
again very poor for J (-0.41) and the P-value very large, 0.071 
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Table 1. Basic properties of the galaxies. 



Galaxy 


R f ,(J) 


R f , (K s ) 


M B 


log V m 


T 


(J) 


y" D (K s ) 


h(J) 


h(KJ 




(Kpc) 


(Kpc) 








(mag/arcsec 2 ) 


(mag/arcsec 2 ) 


(Kpc) 


(Kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 522 


13.7 


11.2 


-20.53 


2.252 


4.1 


20.6 


21.0 


6.6 


7.1 


MGC-01-05-047 


29.3 


27.9 


-21.77 


2.410 


4.9 


20.9 


19.6 


10.8 


8.9 


NGC 781 


11.0 


9.8 


-20.87 




2.4 


19.6 


18.1 


2.9 


2.4 


NGC 2862 


21.7 


18.6 


-21.44 


2.464 


4.0 


19.5 


18.5 


6.0 


5.8 


NGC 3279 


7.8 


6.9 


-19.27 


2.208 


6.8 


19.8 


18.4 


3.9 


3.4 


NGC 3501 


8.2 


7.0 


-19.05 


2.147 


5.9 


20.0 


18.1 


1.9 


1.6 


NGC 4013 


8.5 


7.8 


-19.50 


2.276 


2.9 


18.6 


17.6 


2.2 


1.8 


NGC 4217 


11.1 


10.0 


-19.97 


2.298 


2.8 


18.8 


17.4 


3.2 


2.5 


NGC 5981 


14.6 


12.6 


-20.61 


2.424 


4.7 


19.3 


18.0 


4.2 


3.9 


NGC 6504 


22.2 


23.0 


-22.53 




5.0 


19.5 


18.7 


11.9 


10.3 



Columns: (1) Galaxy name; (2) Truncation radius in J; (3) Truncation radius in K s ; (4) Absolute B-magnitude; (5) Log of maximum velocity 
rotation (in km/s) from radio observations; (6) Morphological type code; (7) Central surface brightness in J; (8) Central surface brightness in 
K s ; (9) & (10) Radial scale-lengths in J and K t . Columns 4, 5 and 6 are obtained from the LEDA database lhttp://leda.univ-lyonl.frl. The 
values of R, r (J), R, r (K s ), h(J) and h(K s ) were independently obtained for both sides of each galaxy. The figure given in this table is the mean 
value. 




h (Kpc) 

Fig. 4. Ratio of truncated radius to scalelength. Symbols as in 
Fig. 1. 



(> 0.05). However, for K s the P-value is less than 0.04 and the 
correlation coefficient -0.73. Therefore, for the K 5 filter there 
is a more significant relation: the larger the y. value, the lower 
the truncation radius when measured with h{K s ) as the unit. 

4. Discussion 

Despite the relatively small size of our sample (10 galaxies), 
some clear relationships between the truncation radius and 
structural parameters of galaxies have been found. These re- 
lations are obtained in the near infrared and therefore some 
caution must be taken when comparing with optical data, as 



Fig. 5. Left panel: Ratio of truncated radius to scalelength ver- 
sus deprojected central surface brightness. Right panel: Always 
considering the scalelength in K s . Symbols as in Fig. 1. 



they could differ. For example, star formation can take place 
beyond the truncation; optical measurements would be sensi- 
tive to newborn stars but not NIR measurements. The difference 
could be very large when comparing data in the U-filter and our 
NIR data and less so when using the R-band and the I-band. It 
should also be taken into account that young stars are visible in 
the NIR, especially red supergiants. Nevertheless, optical and 
NIR data provide different information. This difference would 
be larger if stars moved from their birthplace. The results ob- 
tained are summarized here: 

- Late-type galaxies seem to be more affected by the mecha- 
nism of truncations. 

- Bright galaxies have larger R, r . 

- R, r = constant X V c m (c = 1.5). 

- Larger galaxies (higher h) have larger R, r . 

- < R, r {J)lh{K s ) >= a(J) « 3.57; < R lr (Ks)/h(Ks) >= 
a{Ks) » 3.21; hence, we find lower values for K s . The 
difference is small but R tr (J) > R tr (K s ) for all galaxies with 
only one exception. Statistical tests show that this differ- 
ence is significant. 
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- Rtr/h is an increasing function of the central surface bright- 
ness, I a , even if the correlation seems to be very weak. 

All these properties can be used to put restrictions on the 
different theoretical models. Other groups should confirm these 
results with larger samples of galaxies, with independent NIR 
observations. Clearly, taking into account the fact that some 
galaxies lie behind many Milky Way stars and/or that the ob- 
servations were noisy, the number of galaxies in our sample is 
low and should be enlarged. 

The above results agree very well with the magnetic model 
of Battaner et al. (2002). In this model an inward magnetic 
force is added to gravity to explain the high rotation velocity 
and flat rotation curve of the gas. 

From the observational point of view, it is now evident that 
large magnetic centripetal forces exist, as can be directly de- 
duced from observational magnetic field strengths and energy 
densities in the Milky Way, NGC 6946 and M31 (Beck 2004, 
Battaner & Florido 2005). These measurements show large val- 
ues of the magnetic field strength and, what is more important, 
such low gradients of the field, i.e. such a low decline of the 
field for increasing radii, that an important centripetal magnetic 
force does exist; this is now empirically demonstrated. These 
magnetic forces act on the gas and are in equilibrium with 
gravitational and centrifugal forces. When stars are formed this 
equilibrium is destroyed and stars escape, thus producing the 
truncation. 

In the simplest version of the magnetic model, Battaner et 
al. (2002) obtained a formula for the truncation radius 



2GM 



(1) 



where M is the visible galactic mass and 0q is the constant value 
of the rotation velocity at large radii. 

As M oc L and the Tully-Fisher relation establishes that 
LocV^ (where c lies between 3 and 4), it is concluded that 



(2) 



with the exponent (c - 2) being between 1 and 2. This fact 
matches Fig. 2 very well, showing R, r oc V£, where x = 
1.54 + 0.28 (approximately 1.5 within the estimated error) and 
the constant of proportionality is 0.004 kpc kmT 1 ^ 2 s 1 ' 2 for J. If 
we fit the data in K s the values obtained are: x = 1.58 + 0.29 
and the constant of proportionality is 0.002 kpc kirT 1/2 s 1/2 . On 
the other hand, this would imply that more massive (or more 
luminous) galaxies would have larger R, r , so that the trunca- 
tion might become undetectable if R lr is larger than the last ob- 
servable point. Using the Tully-Fisher relation again, we find 
R tr oc L , c , the exponent (1 - ^) being between 1/3 and 1/2, 
which also explains Fig. 3. 

The magnetic model can also explain the relation between 
R, r /h and /j.o, the surface brightness at the centre. From equa- 
tion (2) and M oc L oc I h 2 (where L is the luminosity and Iq the 
intensity, in physical units, at the centre) and the Tully-Fisher 
relation, taking now for simplicity L oc 



L 1 ' 2 I l l 2 h 



Rtr L 
oc — - — oc oc 

h Vlh h h 



= / l/2 



(3) 



As yU oc - log /, the magnetic model predicts a relation be- 
tween R tr (when measured in radial scale units) and the cen- 
tral surface brightness. Any other characteristic galactic mag- 
nitude (such as h itself) cancels out. 

The magnetic model also predicts R, r /h oc M l/2 /h, and 
hence it explains an observational fact mentioned above: ex- 
tended galaxies, with large values of h, have smaller truncation 
radii (measured with the radial scale length, h, as the unit). In 
this case, the presence of M 1 ^ 2 in the formula would make the 
relation [R tr /h,h] be affected by a high degree of scatter, and 
the relation is certainly more complex than that. 

The fact that the magnetic model satisfactorily explains 
the relations found here does not necessarily imply that all 
other models should be excluded. For instance, on qualitative 
grounds, the threshold model (e.g. Schaye, 2004) could pre- 
dict a relationship between absolute magnitude and truncation 
if the threshold gas surface density were located further out in 
a bright galaxy. Some development of this model to explain 
the findings of this paper is necessary. The magnetic model re- 
mains an interesting alternative, as it predicts not only a quali- 
tative relation, but also the mathematical function. 



5. Conclusions 

The magnetic model satisfactorily explains all the statistical 
properties reported in this paper. We thus propose this scenario 
as responsible for producing truncations. The agreement is par- 
ticularly noticeable concerning the relation between truncation 
radius and the power 3/2 of the asymptotic rotation velocity. 

Magnetic forces are important in the rotation dynamics of 
the gas, as deduced from observed magnetic field strengths. 
These forces no longer act when gas is converted into stars, 
which must have a large influence on the stability of the orbits 
of young stars, even leading to massive escape. 

Large galaxies of type Sa and Sb have less need for dark 
matter, hence, in the magnetic interpretation, they require lower 
magnetic fields, i.e. the truncation would be a less important ef- 
fect in large (massive and luminous) galaxies, especially those 
with a prominent bulge. The fact that late-type galaxies, which 
are bulge-poor, could be more affected by truncations is in 
agreement with the magnetic scenario. 

Pohlen et al. (2004) found several examples of truncations 
in edge-on SO galaxies. Even if these galaxies are at present 
gas-poor, they had both gas and star formation in the past. Thus, 
the mechanisms suggested by the magnetic hypothesis would 
have been operating in the past, and they would have inherited 
truncations. 

The fact that the truncation in the optical could be better de- 
scribed by a two-exponential profile could be compatible with 
a smooth, complete truncation in the near infrared. The newly 
formed stars would be, under the magnetic interpretation, in the 
process of escaping. If the escape speed is roughly equal to V m , 
i.e. of the order of 100 km s _1 , and if the length of the escape 
path is of the order of ~ 10 kpc, the escape time would be of 
the order of 10 8 years, a time that differentiates clearly between 
new-born and old stars. The fact that R t r(K s ) < Rtr( J) also fits 
this interpretation. 
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In the magnetic interpretation, the truncation corresponds 
to a true decrease in the total mass (stars + gas) and a con- 
nection with a decrease in the rotation velocity should be ex- 
pected. This decrease has been observed (Casertano 1983, van 
der Kruit 2001, Bottema 1995) although this fact requires more 
observational support and may be explained by other means. 
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